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Abstract. With the help of 3D MHD simulations we investigate the collapse and fragmen- 
tation of rotating turbulent prestellar core embedded into turbulent medium. The numerical 
code is based on a high resolution Godunov-type finite-difference scheme. Initial turbulence is 
represented by the ensemble of Alfven waves with power law spectrum. Our computations show 
that under realistic parameters two bound fragments can appear when the density increases 
at 100 — 10'^ times. The distance between the fragments is about 0.1 of the initial core radius 
and their orbital period is comparable to the initial free fall time of the core. These results can 
explain the origin of binary stars with separation 0.001 — 0.01 pc in the Galaxy field. 
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1. Introduction 

Prestellar cores in molecular clouds are gravitationally bound condensations with stel- 
lar masses. The collapse and fragmentation of such cores lead to formation of single or 
multiple stars. We investigate the influence of magnetic field, rotation and turbulence 
on a prestellar core collapse and early fragmentation using 3D magnetohydrodynamical 
(MHD) simulations. Similar problem without turbulence was studied by Machida et al. 
(2008), Hennebelle & Teyssier (2008) and Duffin & Pudritz (2009). The role of turbu- 
lence was investigated by Price & Bate (2008) and Wang et al. (2010) for massive clouds 
without rotation. 

We carry out low resolution simulations (with cell number 128'^) and focus the attention 
on initial conditions and results analysis. The computational domain contains both the 
core and surrounding medium to describe a turbulent energy redistribution via MHD 
waves propagation (see Heitsch et al. 2001). The initial pulsations amplitudes depend 
on density (see McKee & Zweibel 1995) because a turbulence in non-uniform medium 
is inhomogeneous one. A bound fragment is found with the help of internal motion 
separation on inward, outward and tangential components relative to fragment mass 
center. 



2. Model 

Initial conditions . Uniform spherical core with radius R — 0.3 of domain size is embed- 
ded into uniform cubic medium (see fig. [l]) with density p — 1/4 of core density. Thermal 
pressure is balanced. Large scale magnetic field Bq is uniform. The core and medium can 
rotate as solid body with axis parallel to Bq. Turbulence is the ensemble of plane Alfven 
waves with linear polarization and isotropic power law spectrum. The initial velocity 
pulsation v-^ in the Alfven wave with wave vector k is directed perpendicularly to (fc, Bq) 
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plane and has amplitude 

Vf:{r) = Afc-"p^(f) sin (kr + (^g) , (2.1) 

where A is the normalization factor, (p^ is the random phase. The magnetic pulsation is 
Bj^ = =FsWj:y/p, where s = signikBo). 

Boundary conditions. Periodical boundary conditions are used to promote the MHD 
waves circulation through the medium. 

Equations and methods. Ideal MHD approximation with isothermal equation of state 
is adopted. The simulations are carried out with the help of numerical code Megalion 
(Dudorov et al. 2003) that realizes the modified Lax-Friedrichs method for hyperbolic 
equations and conjugated gradients method for gravitation equation. The computations 
are single threading. 

Bound fragments. A fragment is the set of contacted cells where density exceeds some 
level. The fragment is bound if 

Eg + Ein > Eth + Em + Etan + Eout , (2.2) 

where Eg is the absolute value of gravitational energy, Etu, Em are the thermal and 
magnetic energies, Eout, Etan are the energies of inward, outward and tangential 
motions relative to fragment mass center. 
Basic parameters. 

1) The ratios of thermal Et, large scale magnetic Em, turbulent kinetic Ek and 
rotational i?n energies of the core to the absolute value of its gravitational energy Eq. 

2) The index of wave spectrum: a £ (1/4, 1/2). 

3) The index of - p relation: /3 e (-1/2,1/4). Note that ^ = -1/2 means the 
balance of wave pressure. 

4) Maximal and minimal dimensionless wave numbers: Umin S (2, 3) and Umax G 
(8,12). 

In all the models Et/Eq = 0.1. The model with parameters Ei^/Eo — 0.05, Ek/Eq = 
0.3, Em/Eg — 0.3, a = 1/3, /3 = —1/2, Umin — 2 and n^ax = 8 is called "base model". 

3. Results 

We have computed successfully 40 models. The basic results are following. 

(a) Turbulence produces density pulsations almost during all the collapse. For example 
at fig. [2] we can see clumpy structure both at the initial and final simulation stages. 

(6) Turbulence usually increases the collapse time by 5 — 10% (fig. E]). Larger (3 leads 
to larger collapse time because of increasing of turbulent pressure gradient. The collapse 
with fragmentation is more delayed than the collapse without fragmentation. 

(c) During the isothermal collapse stage Jeans mass decreases and some fragments 
become bound. Only two bound fragments can appear when the density increases at 
100- 1000 times (fig. Hi). 

(d) The magnetic field directions in the fragments differ from B^ direction. For example 
in fig. |4}d we can see the declinations of magnetic field at 30 — 45° from z-axis. 

(e) In the absence of rotation the early fragmentation is possible when Ek > Et and 
Ek ^ Em (fig. E] left). At moderate rotation {Eq/Eq — 0.05) the early fragmentation 
becomes possible even when Ek = Et and Ek < Em (fig- E] middle) . At faster rotation 
{Eq/Eg = 0.1) the early fragmentation occurs in most cases (fig. [S] right). 

(/) In average, the initial total mass of the bound fragments « 0.07 of core mass, the 
major to minor mass ratio ~ 1.5 (fig. [6]). 
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{g) The distance between the mass centres of the fragments is d — 0.13 — 0.29i?, their 
orbital period is P = 0.47 — 3.2<// (fig.[2l). The correlation between distance and period 
resembles the third Kepler law (P oc cfi^^). 

(h) The properties of bound cores almost independent on initial spectral slope a. This 
results is surprising and must be checked. 

(i) Since typically R ~ 0.01 — 0.1 pc then d ~ 0.001 — 0.01 pc that corresponds to the 
separation in wide binary protostars and stars of the Galaxy field (see Connelley et al. 
2008). 

4. Conclusions 

We carried out 3D MHD simulations of the collapse and fragmentation of rotating 
turbulent prestellar core embedded into turbulent medium. We conclude that: 

• under realistic parameters only two bound fragments can appear when the core 
density increases at 100 — 1000 times, 

• the distance between such fragments is about 0.1 of the initial core radius, 

• their orbital period is comparable to the initial free fall time of the core, 

• these results can explain the origin of binary stars in the Galaxy field with semi- 
major axis in the range 0.001 — 0.01 pc. 

At the next stage of the project we will use the AMR technology realized in the code 
Megalion. That will improve the resolution and lead to more interesting results. It is also 
desirable to develop the turbulent initial conditions without MHD discontinuities at the 
core surface. 
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Figure 1. Initial distributions of density (color scale), velocity and magnetic field (arrows) at 

different planes in the base model. 
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Figure 2. Distributions of density (color scale) and velocity (arrows) at the plane z = in the 
base model at t = 0.144t// and t = 0.964f//. 




Figure 3. Collapse time versus Ek /Eg ratio in all 40 models. Open circles - no fragmentation, 
filled circles - yes. Solid line emphasizes that the collapse with fragmentation is more delayed 
than the collapse without fragmentation. 
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Figure 4. Final distributions of density and velocity at the plane 2: = (a) as well as density 
and magnetic field at the plane y — (b) in the central region of base model. Only density p > 5 
is shown (color scale and contours). 
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Figure 5. Fragmentation diagram: open circles - no, filled circles - yes. See yes-to-no ratio at 

the foot. 
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Figure 6. Fragments mass ratio versus their total mass to core mass ratio. 
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Figure 7. Orbital period of fragments (in core free fall time units) versus their separation to core 
radius ratio. The solid curve shows that the correlation between distance and period resembles 
the third Kepler law (P oc cP^'^). 



